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Abstract 

We show that a recent apphcation of homotopy perturbation method to a class of 
ordinary differential equations yields either useless or wrong results. 

There has recently been great interest in the application of several approximate 
procedures, like the homotopy perturbation method (HPM), the Adomian 
decomposition method (ADM), and the variation iteration method (VIM), to 
a variety of linear and nonlinear problems of interest in theoretical physics 
[1-15]. For brevity I will call VAPA all those variational and perturbational 
approaches. In a series of papers I have shown that most of the results produced 
by those methods are useless, nonsensical, and worthless [16-19]. 

In a recent paper Rafiq et al [13] applied the HPM to some ordinary second- 
order differential equations, and the purpose of this comment is to discuss 
their results. 
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Rafiq et al [13] solved second order differential equations of the form 

y"{x) +p{x)y'{x) + f{x,y) = 

(1) 

y{0) = A, y'{0) = B 

where /(x, y) = F{x, y) — g{x) in their notation. In order to apply HPM they 
wrote y"{x) +p{x)y'{x) + 9f{x, y) = and expanded the solution in a 6'-power 
series 

y{x) = yo{x) + Oyiix) + . . . (2) 

Finally, they set ^ = 1 in order to obtain an approximate solution to equation 
©• 

In the particular HPM implementation proposed by Rafiq et al [13] the per- 
turbation corrections yj{x) result to be polynomials so that the partial sums 
of the HPM series are just x-power series of the form 

y{x) = A + Bx + a2X^ + a^x^ + ... (3) 



In what follows we analyze the examples chosen by Rafiq et al [13]. The first 
one is a textbook exercise for an introductory course on differential equations: 

y"{x) + -y'{x) +xy = x^ + AAx"^ - 30x 

(4) 

A = B = 

We appreciate that y{x) ~ as x -^■ oo. Therefore, if we substitute the 
polynomial y{x) = a2X^ + a^x^ + a^x^^ into equation (jl]) we easily obtain the 
exact solution y{x) = x'^ — x^. However, Rafiq et al [13] apply a cumbersome 
perturbation method and obtain an infinite number of nonzero perturbation 
corrections yj{x) that cancel out to (hopefully) produce the exact solution. 
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The authors do not prove that their HPM aheady yields the exact result, 
and any partial sum gives the exact result plus terms that are cancelled by 
corrections of higher order. In this case the HPM partial sum of any order is 
an inexact approach to an extremely simple solution of a trivial differential 
equation that one easily derives by the straightforward method just indicated. 

The second example is as trivial as the first one: 

y"{x) + -y'{x) + y = 6 + 12a; + 

(5) 

A = S = 

Any undergraduate student will try y{x) = a2x'^ + a^x^ and obtain the exact 
solution y{x) — x'^+x^ without effort. Again, the HPM yields an infinite series 
and the authors do not prove their convergence. Although it seems that the 
spurious terms cancel out, any partial sum yields a wrong result. 

A simple inspection of the third example 
y"{x) + + = 6 + 

(6) 

suggests that y{x) — a2X^ and thus one obtains the exact result y{x) — x'^.By 
means of the HPM Rafiq et al [13] obtain the wrong result y{x) = x'^ + x^/72\\\ 
Chowdhury and Hashim [4] also obtain this wrong result but then they choose 
a different starting point of the perturbation approach in order to derive the 
expected exact solution. In the first two examples Rafiq et al [13] mention the 
appearance of noise terms. It seems to me that it was that noise that already 
affected the calculation in this example. 
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The fourth example is much more interesting (at least it is not trivial): 
y"(x) + -y'ix) + e^?^' = x + 1 

(7) 

A = S = 

After a long and tedious perturbation calculation the authors derive the first 
correct terms of the Taylor series expansion about a; = [13] 

y(x) = - ^ + . . . (8) 

^ 12 12960 ^ ^ 

It is unnecessary to say that one can easily obtain the same approximate 
result more easily by means of the power-series method. We expect that the 
power series will not reveal the most interesting features of the solution to this 
equation but a mere indication of what happens in a neighbourhood of a; = 0. 

At first sight one guesses that the solution to equation ([7]) should behave 
approximately as y{x) ~ ^\ia{x)/x as x ^ oo; however, its actual behaviour is 
richer. Fig. [1] shows that y{x) (calculated numerically with sufficient accuracy) 
oscillates about ^\n(x)/x and approaches this asymptotic function as x ^ oo. 
The x-power series ([8]) only accounts for the behaviour of y{x) up to about 
the first maximum. If equation ([7]) represented an actual physical problem we 
would be missing its most interesting features when using the HPM to obtain 
its solution. 

Summarizing: the HPM proposed by Rafiq et al [13] gives cumbersome approx- 
imations to simple solutions of trivial differential equations, yields a wrong 
result in one of the cases, and fails to provide the most interesting features 
of the solution to the only nontrivial example. Unfortunately the authors do 
not indicate any physical application of the examples chosen. They seem to 
be just toy problems for fiddling around with the HPM. 

The conclusion above is not surprising because all our previous analysis of 
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several applications of VAPA have shown that they produce useless or trivial 
results [16-19]. See, for example, the earlier papers by Chowdhury and Hashim 
[4] and Yildirim and Ozi§ [3] who applied similar cumbersome perturbation 
equations and obtained the power series for the solutions of exactly solvable 
differential equations. 

It seems that one of the greatest feats of many VAPA applications is to produce 
power-series expansions to simple and trivial problems in a cumbersome and 
laborious way [3,4,13]; the reader may find the analysis of other such examples 
in our earher communications [16-19]. In fact, VAPA have produced the worst 
research papers ever written. 
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Fig. 1. Exact (accurate numerical) solution, power series ([8]) and asymptotic expan- 
sion for example ([7]) 
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